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Abstract 



We have generalized the network approach to include the effects of short- 
range imperfections in order to analyze recent experiments on mesoscopic 
superconducting double loops. The presence of weakly scattering imperfec- 
tions causes gaps in the phase boundary B(T) or <£(T) for certain intervals of 
T, which depend on the magnetic flux penetrating each loop. This is accom- 
panied by a critical temperature T c ($), showing a smooth transition between 
symmetric and antisymmetric states. When the scattering strength of imper- 
fections increases beyond a certain limit, gaps in the phase boundary T C (B) or 
T c (<3?) appear for values of magnetic flux lying in intervals around half-integer 
$0 = hc/2e. The critical temperature corresponding to these values of mag- 
netic flux is determined mainly by imperfections in the central branch. The 
calculated phase boundary is in good agreement with experiment. 
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Early experimental! have revealed that the effect of nonmagnetic impurities on the 
transition temperatures of bulk superconductors is very small. The critical temperature, 
T c , changes by about 1% for 1% concentration of impurities. An interpretation of these 
observations was first proposed by Anderson!. The only effect of the impurities is to change 
the energy of a free electron to eigenenergies determined by those impurities. This modifies 
the density of states in the integral equation for T c . Therefore, scattering by non-magnetic 
impurities only slightly changes T c in bulk superconductorsi. Investigations extending An- 
derson's work have been performed for multiband superconductors using the Abrikosov- 
Gor'kov approachi. An interesting opportunity to intensify the effect of the imperfections 
occurs in mesoscopic superconducting structures where the confined condensate is much 
more sensitive to the action of impurities than in bulk structures. 

Recently, the onset of the superconducting state has been studiedl'0 in different meso- 
scopic structures of Al, comprising lines, dots, loops, double loops, microladders etc., with 
sizes smaller than the coherence length £(T). Refs. |S]|| show that experimentally observed 
phase boundaries for square loops with two attached leads are in excellent agreement with 
calculations based on Ginzburg-Landau (GL) theory0'0. 

In this present letter, we analyze the case of a superconducting In this present letter, 
we analyze the case of a superconducting mesoscopic double square loop (see the inset to 
Fig. la) where experiments § have revealed the phase boundary shown in Fig. la. Using 
the micronet approachEHlil for a superconducting double loop, we obtain a phase boundary 
(cf. Ref. |7|) which consists of the intersecting parabolas shown in Fig. lb. One set (with 
minima at integer values of the magnetic flux $ through one loop in units of the magnetic 
flux quantum $ = hc/2e: $/$o) depends on the magnetic flux quanta penetrating each 
loop with L = 0,1,2,.... The other set (with minima at half-integer values of $/$o) 
depends on odd numbers of magnetic flux quanta penetrating the double loop as a whole 
with L — 1/2, 3/2, . . .. Since the critical temperature corresponds to the lowest Landau level 
Elll{&), the way the parabolas intersect means that the minimum energy encounters a shift 
from one branch to another at certain values of magnetic flux. Moreover, since the derivative 
of the lowest Landau level with respect to magnetic flux is proportional to the persistent 
current (cf. Ref. [15| Eq. (4.5)) we arrive at the following paradox: At the intersections, the 
left and right derivatives of Emi^) are different, the persistent current has a discontinuity 
and its value is consequently indefinite. In order to resolve this paradox, an analogy between 
superconducting loops and semiconducting quantum rings can be exploited. In such rings, 
in the presence of impurities and an Aharonov-Bohm magnetic field, a crossing is known to 
change into an anti-crossing. In this case, the gaps between the different eigenenergies as a 
function of the magnetic flux widen and hence the degeneracy of states, which contribute to 
the lowest level, is raised (see Ref. [16]). In this letter we demonstrate that this also applies 
to the superconducting mesoscopic double square loop. 
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FIG. 1. Experimental phase boundaries for a double loop of Al (shown in the inset) with (panel 
a) and without (panel b) a parabolic background (the dashed line in panel a) due to the finite width 
of the loop. When transforming the phase boundary to the plane T c ($)/T c versus &/&o, the value 
Q=1302 nm is taken in order to ensure that the maxima of T c (&)/T c are at integer values of ^/^o- 
Theoretical phase boundaries (solid lines in panel b) calculated for a perfect double loop with 
Q=1302 nm and £ = 128 nm. 

The observed phase boundary for a superconducting mesoscopic double square Al loop is 
given in Fig. la. After subtraction of the parabolic background, related to the finite width, 
w = 130 nm, of the stripes (dashed line in Fig. la), at least 12 practically identical periods 
are seen. One such period is plotted in Fig. lb (dots). From the period of the oscillations 
of the phase boundary with respect to the magnetic field, AB = 1.24 mT, an effective loop 
side length Q = 1.3 /iin is obtained. This is close to the average loop size. For further 
experimental details, the reader is referred to Ref. []. In order to understand the observed 
"anti-crossing" of the different elements of the experimentally observed phase boundary 
and, in particular, the smooth shape of the T c ($) minima, we consider a superconducting 
mesoscopic double square loop with imperfections. These imperfections may be introduced 
during the fabrication of the mesoscopic structures. Most probably one of the sources of 
such imperfections is the inhomogeneity of the geometrical superconducting lines written by 
e-beam lithography. 

Within the framework of the GL approach, the presence of imperfections in a supercon- 
ducting structure may be modelled by a spatial inhomogeneity in the parameters a and b in 
the GL equation: 

— (-iKV - — A(r)) V(r) + a(rWr) + 6(r) I ^(r) I 2 ip(r) = 0. (1) 
2m V c / 

Near the phase boundary, where the order parameter ip(r) is small, the system is adequately 
described by the linearized GL equation: 

— f-inV - — A(r)) Mr) + a(rU(r) = 0. (2) 
2m \ c J 

Moreover, the magnetic field may be assumed to be equal to the applied magnetic field. The 
vector potential A of the uniform field B || e z is taken in the symmetric gauge. 

The presence of imperfections, localized in the loop around several points r s over a 
distance which is much smaller than the coherence length or the typical loop size ("short 
range imperfections"), can be modelled by the following function: 
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a(r) = a + Y,V s 5(r-r s ), (3) 

s 

where a is the GL coefficient of the substance, and the magnitudes V s are determined by 
specific characteristics of the imperfections. Eq. (0) then becomes similar to the Schrodinger 
equation for a particle of mass m and charge 2e in the potential field described by the scalar 
form J2 S K<5(r — r s ) and by the vector potential A(r), the quantity —a playing the role of 
the energy. 

Short-range imperfections are assumed to be present in all three branches of the loop 
at the points characterized by the coordinates Q s with s = L, M, R as shown in Fig. 2. 
Furthermore, taking £(T) = £o/\/(l — T/T c ) as a unit of length, we obtain the linearized 
GL equation in terms of dimensionless coordinates: 

( 2ir 2tt - 1 

(, Vx + —A x f + (<V y + —A y f + £ K%-Q,)-lk = 0. (4) 

I ^0 S =L,M,R > 

It should be noted that the dimensionless scattering magnitude V s = 2mV s ^(T) /h 2 = 
C 8 / J(l — T/T c ) is temperature dependent. 
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FIG. 2. A model configuration of imperfections (open circles) in the double loop. 



We have solved Eq. (f|) following the micronet approach. A new conceptual feature 
compared with Refs. [12|-|14| is the use of additional nodal points at the positions of the 
imperfections. The presence of these imperfections implies an additional condition, which 
can be derived in the following way. In the vicinity of the point with y = Q s , Eq. (|J) takes 
the form 



2-7T 

(tV y + ^-A y ) 2 ^ - * = -V s 5(y - Q a )tt. 

We integrate both sides of this equation over y from Q s — e to Q s 
account continuity of the order parameter, ^(Q s + e) = ^(Q s — e) 
constraint 



(5) 



■f e, e — > +0. Taking into 
we obtain the additional 



dm 

dy 



y=Q s +e 



dm 

dy 



V s m(y = Q s 



(6) 



y=Qs 



for the derivatives to the left and to the right from the new nodal point y = Q s . Applying 
Eq. (||) to a single loop of length C with one imperfection leads to the following secular 
equation 
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cos((p) = cos(£) + -V s sin(£), 



(7) 



where tp = 27r$/$ - 

Analyzing the onset of superconductivity in a superconducting ring with a lateral arm 
of length £, de Gennest3 derived an equation which differs from Eq. (|7|) by the substitution 
V s = — tan(£)/2. Thus, from the point of view of our present analysis, a lateral arm may be 
considered of CIS db kind of imperfection in a ring. 

After some lengthy, but straightforward algebra, we obtain the secular equation which 
determines the phase boundary: 
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with 

N u = cos3Q + V s sin(2Q - Q s ) cos(Q + Q s ), N 2s = cos3Q + V" s cos(2(5 - Q s ) sin(Q + Q s ), 
D s = sin3Q + V s sin(2Q - Q s )sin(Q + Q s ); 

Nim = cosQ + Vm sin(Q - Qm) cos(Qm), N 2 m = cosQ + V M cos(Q - Qm) sin(Q + Q M ), 
D M = sinQ + V M sm(2Q - Q M ) sin(g + Qm), (9) 

where s = L,R. 

In Eq. (^[), ip — 27T$/$ with $, the magnetic flux through each of the loops. Here, we 
recall that lengths in the above equations are expressed in units of £(T) and are therefore 
functions of the temperature. Consequently, the secular equation establishes a relation 
between the magnetic flux $ and the temperature T. 

The phase boundaries obtained by solving Eq. (f|), with an imperfection in only one 
branch of the double loop, are shown in Figs. 3 to 5 {L in Figs. 3 and 4, M in Fig. 5). 
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FIG. 3. Theoretical phase boundaries for a double loop with an imperfection in one side branch 
at Ql = 0.5Q for (5=1302 nm and £o = 128 nm. (a) The dimensionless magnitude Cl= and 
0.01 corresponds to dotted and solid lines, respectively, (b) The dimensionless magnitude Cl= 
0, 0.5, 1.0 corresponds to dotted, solid and dashed lines. The values $i and $2 are discussed in the 
text. 
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As is seen in Fig. 3a, at small values C s , an "energy" gap forms between solu- 
tions corresponding to different numbers of magnetic flux quanta penetrating each loop 
(L = 0, 1/2, 1, ...). The resultant phase boundary indicates that a continuous change takes 
place from a symmetric superconducting order parameter at integer values of $/$o to an 
antisymmetric state at half-integer values of the relative magnetic flux. It is also worthy of 
note that the presence of imperfections slightly diminishes the critical temperature of the 
double loop at zero magnetic field. 

When increasing C s above a certain limit, the pattern of phase boundaries changes 
dramatically. Gaps appear in certain flux intervals ("flux gaps") around half- integer $/$o 
values. This behavior is illustrated by the curves in Fig. 3b. For a given T < T c , a 
superconducting state exists when $ ranges from zero to a certain value of $i, after which 
the sample turns into the normal state. With a further increase of magnetic flux along a 
horizontal straight line shown in Fig. 3b, the sample remains in the normal state until a value 
of $2 is reached, at which the sample becomes again superconducting. This demonstrates 
a re-entrant behavior as a function of field. [When approaching in Fig. 3b the points where 
the phase boundary <£>(T) has its extrema and, consequently, the derivatives <9T/cM> would 
diverge, the superconducting state apparently becomes unstable.] It should be noted that 
the existence of such a regime, where the system is normal in a certain flux interval, was 
reported by de Genneslll for a single superconducting ring with a lateral arm. 

The trend of lowering of the critical temperature of a double loop at zero magnetic field 
with increasing C s is clearly seen by comparison of the curves in Fig. 3b which refer to 
different values of Vl- 




FIG. 4. Theoretical phase boundaries (solid lines) for a double loop with an imperfection in 
one side branch at Ql = 0.5(3 f° r Q=1302 nm and £o = 128 nm. The dimensionless magnitude is 
Cl= 0.02. Results obtained without and with renormalization of temperature are given in panels 
a and b, respectively. Experimental data are shown with points. 

In regions adjacent to integer values of $/$cb a good agreement with experiment of the 
calculated phase boundary T c ($), is achieved for the zero-temperature coherence length £ 
= 128 nm. This value is in accordance with previous estimates (see Ref. [17]). The resulting 
set of phase boundaries is shown in Fig. 4a. For comparison with the experimental data, it 
is necessary to renormalize the temperature scale, taking as a unit temperature the specific 
critical value of a loop with imperfections (see Fig. 4b). 

In Fig. 5, a plot of the phase boundary is shown for a double loop with an imperfection 
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in the middle branch (M in Fig. 2). It is clear that by increasing Vm, one shifts a minimum 
of the 1 — T c ($)/T c curve (using a renormalized temperature), at half-integer values of $/$o? 
to higher values, without modifying those parts of the phase boundary close to the integer 
values of $/$ - 
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FIG. 5. Theoretical phase boundaries for a double loop with an imperfection in the middle 
branch at Qm = 0.1Q for (5=1302 nm and £o = 128 nm. The dimensionless magnitude Cm= 
0.11,0.125 corresponds to dashed and solid lines. Results are obtained with renormalization of 
temperature. Experimental data are shown with points. 

The best agreement between the calculated phase boundary and experimental data is 
achieved for a configuration where imperfections are present on all three branches of the 
double loop (L,M,R in Fig. 2). The corresponding curves are shown in Fig. 6. The main 
observation which follows from these figures is that imperfections in the central branch play 
a decisive role in determining the critical temperature of a double loop at the half-integer 
values of $/$ - 
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FIG. 6. Theoretical phase boundaries (solid line) for a double loop with imperfections in all 
three branches at Ql = Qr = Qm = 0.5Q for (5=1302 nm and £o = 128 nm. The dimensionless 
magnitudes are: Cl= —Cr = 0.017 and Cm= 0.090. Results are obtained with renormalization 
of temperature. Experimental data are shown with points. 



In conclusion, we have generalized the network approach to include the effects of short- 
range imperfections. The presence of weakly scattering imperfections leads to the formation 
of gaps between solutions corresponding to integer and half-integer numbers of magnetic 
flux quanta penetrating each loop. The phase boundary is characterized by a smooth tran- 
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sition between symmetric and antisymmetric states. For imperfections with relatively large 
magnitudes, gaps in the phase boundary T C (B) or T c (f ) appear when the magnetic flux lies 
in intervals around half-integer values. The critical temperature at the half-integer values 
of the relative magnetic flux has been shown to be determined mainly by imperfections in 
the central branch. The calculated phase boundary for a mesoscopic double square loop is 
in good agreement with experiment. 
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